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The Pierre Auger Observatory's (PAO) shower profile measurements can be used to constrain the 
chemical composition of the ultra-high energy cosmic ray (UHECR) spectrum. In particular, the 
PAO's measurements of the average depth of shower maximum and the fluctuations of the depth of 
shower maximum indicate that the cosmic ray spectrum is dominated by a fairly narrow distribution 
(in charge) of heavy or intermediate mass nuclei at the highest measured energies {E ^ 10^^ eV), and 
contains mostly lighter nuclei or protons at lower energies {E ~ 10^* eV). In this article, we study 
the propagation of UHECR nuclei with the goal of using these measurements, along with those 
of the shape of the spectrum, to constrain the chemical composition of the particles accelerated 
by the sources of the UHECRs. We find that with modest intergalactic magnetic fields, 0.3 nG 
in strength with 1 Mpc coherent lengths, good fits to the combined PAO data can be found for 
the case in which the sources accelerate primarily intermediate mass nuclei (such as nitrogen or 
silicon). Without intergalactic magnetic fields, we do not find any composition scenarios that can 
accommodate the PAO data. For a spectrum dominated by heavy or intermediate mass nuclei, the 
Galactic (and intergalactic) magnetic fields are expected to erase any significant angular correlation 
between the sources and arrival directions of UHECRs. 



I. INTRODUCTION 

The chemical composition of the ultra-high energy cos- 
mic ray (UHECR) spectrum has long been a topic of 
great interest [THS]- Until recently, however, very little 
was known about the nature of these particles. On one 
side of the debate, the so called Hillas criterion T gives 
a preference for the electromagnetic fields of cosmic ray 
sources to accelerate heavy nuclei to higher energies than 
protons or light nuclei. On the other side, it has been ar- 
gued that the angular correlations reported by the Pierre 
Auger Observatory (PAO) [5], as well as features in the 
shape of the UHECR spectrum [5], suggest that these 
particles consist largely of protons. None of these ar- 
guments, however, has yet settled the question of what 
types of particles make up the UHECR spectrum. 

Data from the Pierre Auger Observatory (PAO), how- 
ever, is offering increasingly powerful insights into this 
question. Firstly, the spectral shape predicted for the 
UHECR all-particle spectrum depends not only on the 
injected spectrum and spatial distribution of the sources, 
but also on the chemical composition that is injected 
from the sources of the highest energy cosmic rays. As 
the PAO measures the UHECR spectrum with increasing 
precision ^lOj . this information can be used to constrain 
the chemical composition of these particles [TT] . Further- 
more, the PAO is capable of performing several measure- 
ments that can be used to directly or indirectly deter- 
mine the chemical composition of UHECRs as they enter 
the Earth's atmosphere. Among these empirical tools 
are the measurements of the average depth of shower 
maximum, (Xmax) , and the RMS variation of this quan- 
tity, RMS(A'i„ax)- On average, proton-induced showers 



reach their maximum development, (v'Cmax), deeper in 
the Earth's atmosphere than do showers of the same en- 
ergy generated by heavier nuclei. Accompanying this re- 
sult, the shower to shower fluctuation of A"max about the 
mean, RMS(Xmax), is larger for proton-induced showers 
than for iron-induced showers of the same energy. As 
a result, measurements of both (Xmax) and RMS(Arinax) 
can be used to infer the average chemical composition of 
the UHECRs as a function energy. 

Very recently, the PAO collaboration has announced 
their first measurements of RMS(Xniax) [HI US]- These 
measurements, along with those of (A^max), imply that 
the UHECR spectrum contains a large fraction of heavy 
or intermediate mass nuclei, especially at the highest 
energies measured. Furthermore, the small values of 
RMS(X,nax) measured by the PAO also imply that the 
composition of the UHECR spectrum is relatively nar- 
rowly distributed at the highest measured energies, con- 
taining little or no protons or light nuclei. In this way, 
the new RMS(Xinax) measurements not only confirm 
and reinforce the conclusions drawn from earlier aver- 
age depth of shower maximum measurements, but also 
provide complementary information that enables one to 
constrain the distribution of the various chemical species 
present within the UHECR spectrum. 

The remainder of this article is structured as follows. 
In Sec. [IT] we review the physics of UHECR nuclei prop- 
agation. In Sec. |III[ we calculate the spectrum and 
chemical composition of the UHECR spectrum for var- 
ious choices of the injected chemical species and com- 
pare these results to measurements of the PAO, ne- 
glecting the effects of intergalactic magnetic fields. In 
Sec. |IV[ we include magnetic fields in our calculation, 
and show that fields with stengths and coherent scales 
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'-^(i?/0.3 nG)x(Lcoh/l Mpc)-^/^, are strongly preferred to 
accomodate the PAO's data. The best fits found are for 
cases in which the sources of the UHECRs inject mostly 
intermediate mass nuclei, such as nitrogen or silicon. Fi- 
nally, in Sec.|Vj we discuss the implications of our results 
and summarize our conclusions. 



II. THE PROPAGATION OF ULTRA-HIGH 
ENERGY COSMIC RAY NUCLEI 

For cosmic rays in the form of ultra-high energy (UHE) 
protons, the dominant processes effecting their propa- 
gation over cosmological distances are their interactions 
with the cosmic microwave background (CMB) or cos- 
mic infrared background (GIB), producing either pions 
or electron-positron pairs. Pair production (p -t- 7 
p + e+ -|-e~) dominates for energies below lO^^'® eV, and 
occurs with sufficiently small inelasticity that it can be 
treated as a continuous energy loss process [H] . Pion pro- 
duction, by contrast, occurs with large inelasticity, with 
individual occurrences of the processes p-l-7CMB — >■ P+tt^ 
or p-l-7cMB n-t-Tr"*" causing the primary proton to lose 
a considerable fraction of its energy. Thus pion produc- 
tion must be treated as a stocastic process, and follow- 
ing its onset leads to the catastrophic attenuation of the 
UHECR spectrum above ~ lO^^-^ eV known as the GZK 
cutoff [15\ Furthermore, if enough center-of-mass en- 
ergy is available, multi-pion production can also become 
important. 

The propagation of UHE cosmic ray nuclei is some- 
what more complicated. In addition to energy losses from 
pair-production, cosmic ray nuclei undergo photodisinte- 
gration in scattering off the GMB and/or GIB at a rate 
given by: 

(1) 

where A and Z are the atomic number and charge of the 
nucleus, ip and in are the numbers of protons and neu- 
trons broken off from a nucleus in the interaction, 71(e) 
is the density of background photons of energy e, and 
fA,z,ip,i„(e') is the appropriate cross-section [TH31IT6]. In 
our calculations, we use the infrared background model 
of Ref . J/Zj , although other models [TS] yield similar re- 
sults. As a result of the process of photodisintegration, 
the chemical composition that reaches Earth can be sig- 
nificantly different than the composition that is initially 
injected from the sources of UHECRs. 

In previous work, we have described our Monte Garlo 
code which models the propagation of UHE cosmic ray 
nuclei over cosmological distances (see Ref. [3]). In this 
article, we use this code to determine the spectrum and 
composition of UHEGRs that is predicted to reach Earth, 
and compare this to the measurments of the PAG for a 
variety of injected chemical compositions. 



III. RESULTS IN THE CASE OF NEGLIGIBLE 
MAGNETIC FIELDS 

To begin, we consider the case in which only iron nuclei 
are accelerated by the sources of UHECRs, and assume 
that the deflection of UHEGRs by intergalactic magnetic 
fields is of negligible importance. We adopt a homoge- 
neous distribution of sources, and a simple cosmic ray 
injection spectrum of the form: 

dE ^ ' 

It is important to recognize that this describes the spec- 
trum that is produced by the sources of the UHECRs, 
which can be very different from the spectrum (and com- 
position) that exists after the effects of propagation are 
taken into account. As a function of the spectral slope 
injected, a, and the energy cutoff, -Emax, we can calcu- 
late the resulting spectrum and composition predicted to 
be observed at Earth, and compare this to the measure- 
ments of the PAO. 

In the following we obtain best fits to the arriving spec- 
trum, (Xniax), and RMS(Armax) data by scanning over 
Emax from 10^°-10^^ eV, and a from 1.6-2.4, and min- 
imising the for these runs, keeping the best fit results. 
Only data for GR energies above 10^*'^ eV were used in 
these fits, with the assumption that events below this 
energy are Galactic in origin. All spectral results are 
shown for energies above this assumed transition energy. 
However, to provide a comprehensive understanding of 
the change in composition following photodisintegration 
interactions, all composition related plots ((^max), and 
RMS(Xmax)) are shown down to energies of 10^^ eV. We 
reiterate to the reader that only data points with ener- 
gies above 10^*'^ eV have been used to obtain the best 
fit results. In order to keep the composition related plots 
as simple as possible, we plot only the upper and lower 
(Xmax) and RMS(Xinax) values allowed from the different 
hadronic models. For each shaded region in the (ATmax) 
plot the upper and lower boundaries are predominantly 
dictated by the EPOS and QGSJET models respectively. 
No simple rule, however, explains the model responsible 
for the upper and lower boundaries of the shaded regions 
in the RMS(Ainax) plot. 

In Fig. [1] we show our results for the cases of i?niax — 
lO^O-S and lO^i eV (a =1.6 in each case). We note that 
although such hard spectra are not well motivated by 
non-relativistic diffusive shock acceleration theory, such 
spectra have been suggested by relativistic shock accel- 
eration theory |19j . Furthermore, the reader should bear 
in mind that such hard spectra obtain good agreement 
through the considerable lightening of the arriving com- 
position that they allow, caused through the additional 
contribution of liberated protons following photodisinte- 
gration interactions. Thus hard injection spectra provide 
the possibility of a lightened composition at lower ener- 
gies. This lightening of the composition, however, comes 
at the price of a reduced spectral fit with the PAO data. 
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FIG. 1: The spectrum (left), the average depth of shower maximum (center), and the RMS variation of the depth of shower 
maximum (right), as measured by the Pierre Auger Observatory (PAO) [101 112| . compared to the predictions in the case that 
the sources of the ultra-high energy cosmic rays accelerate only iron nuclei. The iron nuclei were assumed to be injected from 
their sources with a spectrum of dN/dE cx E~°'e~^^^'^'-^ , where -Emax is 10^"'^ or 10^^ eV in the upper and lower frames, 
respectively. In each case, the spectral index a was selected to provide the best fit to the spectrum measured by the PAO 
(a = 1.6 in each case shown). In the center and right frames, we show the regions bounded by the predictions of four different 
hadronic models (QGSJET [22], SIBYLL |23) . and EPOS [24]), and for comparison also show each simulation's prediction for 
a spectrum at Earth consisting of purely protons or purely iron nuclei. For each shaded region in the shower maximum plot 
the upper and lower boundaries are predominantly dictated by the EPOS and QGSJET models respectively. No simple rule, 
however, explains the model responsible for the upper and lower boundaries of the shaded regions in the RMS variation plot. 



Such good fits, for the case of heavy nuclei, therefore 
exist through the balance between these two effects. In 
the middle and right frames, we have used the results 
of hadronic simulations to convert the chemical compo- 
sition of the spectrum found by our Monte Carlo into 
the related quantities observed by the PAO - the aver- 
age depth of shower maximum, and the RMS variation 
of this quantity The shaded region in the plots is 

that bounded by the results obtained using four differ- 
ent commonly used parameterizations |13) of air shower 
simulations [21]: QGSJET II, QGSJET [22], SIBYLL 
2.1 [23], and EPOS [24 . For comparison, we also show 
the results that are obtained for the cases in which ei- 
ther only protons or iron nuclei are present in the cosmic 
ray spectrum at Earth. At some level, the variations 
between the results of the different hadronic simulations 
provide us with an estimate of the systematic errors asso- 
ciated with our ignorance of the non-perturbative QCD 
necessary for a complete description of the shower de- 
velopment. Furthermore, we highlight to the reader the 
fact that these models only provide an incomplete scan 
of the full parameter space [25H28] . leaving room for the 
possibility that the model actually realised by nature sits 
outside this range. 

From Fig. [ij we see that in the case of all-iron in- 



jecting sources with E'max It, lO^"'^ eV, the measured 
average depth of shower maximum consistently exceed 
the predicted values. This is easy to understand from 
the fact that 10^°'^ eV iron nuclei fragment into nucle- 
ons with only ^ 6 x 10^® eV of energy each, thus leaving 
only heavy nuclei to make up the spectrum at higher en- 
ergies. As the maximum energy is increased, however, 
more protons are present among the UHECRs, causing 
the average composition to be considerably lighter. If we 
increase E'max, we can better accomodate the measure- 
ments of (Xmax), but also begin to significantly exceed 
those of RMS(Arnjax)- For an all-iron source composition, 
we do not find any combination of ii'max and a that can 
simultaneously accomodate the PAO's measurements of 
(^max) and RMS(X 

max / ■ 

To see whether less heavy species of cosmic ray nuclei 
might be able to better provide a reasonable fit to the 
PAO data, we show results in Figs. |2] and |3] for the cases 
of a pure silicon or pure nitrogen injection spectrum. For 
the all-silicon case, the best fit corresponds to the choices 
of -E-max ~ 10^^ eV and a 1.8.^ Even in this best 



Throughout this study, the quantity -Bmax refers to the cutoff 
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FIG. 2: As in Fig.[T] but for the case in which only sihcon nuclei are accelerated by the sources of the ultra-high energy cosmic 
rays. 
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FIG. 3: As in Figs. [T] and [2] but for the case in which only nitrogen nuclei are accelerated by the sources of the ultra-high 
energy cosmic rays. 



case, however, the fit to the (ATmax) data is quite poor, 
although better than that found in the all-iron case. 

The all-nitrogen case even more clearly fails to accom- 



energy for iron nuclei. This is related to the cutoff for other 
species of charge Z by Emax.z 

= iZ/26) X 



modate the observations, in particular with regards to 
the observed shape of the UHECR spectrum which, un- 
like the composition measurements, are not impacted by 
uncertainties associated with the hadronic models. This 
is in agreement with the conclusions drawn in Ref. [llj . 
Furthermore, the all-nitrogen case does not predict a suf- 
ficiently rapid change of RMS(Armax) with energy. We 
thus conclude that while the PAO observations do not 
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appear to be very well fit by any of the scenarios we have 
considered so far, the best {ie. least bad) fits are in those 
cases in which the observed UHECR spectrum originates 
largely from sources accelerating intermediate mass or 
heavy nuclei {A ^ 20). 

If cosmic ray sources accelerate a mixture of dif- 
ferent species of nuclei, this will impact (Xmax) and 
RMS(Xiiiax) in somewhat different ways. In particu- 
lar, whereas the average depth of shower maximum at a 
given energy simply scales with (Xj^ax) oc 1/ (In A), the 
RMS(Xmax) for an ensemble of different nuclear species 
is given by, 

56 

O'tot = X! (f^'^A + /A(^max,A - (^max))^) , (3) 
A=l 

where /a is the fraction of the total population with 
atomic mass A, and (ta is the RMS of X^^x predicted for 
cosmic rays of that atomic mass. From this, we see that 
the relatively low value of RMS(Xijiax) observed by the 
PAO not only indicates that a large proportion of the ar- 
riving nuclei are heavy at the highest energies, but also 
that their composition is fairly narrowly distributed in 
charge- ie. the fraction of protons or light nuclei in the 
UHECR spectrum must be small in order to prevent large 
contributions from the /i(X„iax,i ~ (^max))^ terms in the 
sum. In many cases, this leads to a tension between the 
PAO's measurements of (Xmax) and RMS(Xmax), which 
can be made worse by considering cosmic ray sources 
which inject particles with a mixed chemical composition. 
We find that mixed (more than 2 component) composi- 
tion scenarios, suggested to be indicated by the recent 
PAO data [29], do not lead to sufficiently improved fits 
to justify their consideration. 

In Fig. |4j we show results for some selected cases of 
mixed chemical composition. In particular, we consider 
combinations of protons, nitrogen, silicon, or iron, and, 
as before, have selected in each case the injected spec- 
tral index, energy cutoff, and fractional mixture that 
provides the best fit to the PAO spectrum, (Xmax) and 
RMS(Xi„ax) measurements. We find that while either the 
(-''^max) or the RMS(X„iax) data can be accommodated 
without much difficulty, it is again challenging to find an 
injected spectrum and composition that is in agreement 
with both (X,nax) and RMS(Xniax) simultaneously. In 
particular, if (X„iax) is well fit, then the predictions for 



RMS(Xmax) tend to exceed the measured values. Simi- 
larly, if the measurements of RMS(Xi„ax) are well fit, the 
predictions for (X^ax) tend to be well below the mea- 
sured values. The only exception we find to this conclu- 
sion is the case in which the injected spectrum consists 
of a mixture of nitrogen nuclei (80%) and protons (20%), 
which (using the EPOS hadronic simulation) can pro- 
vide a reasonable fit to both (Xmax) and RMS(Xniax) 
measurements. In this case, however, the spectral shape 
drops far more rapidly than measured by the PAO, lead- 
ing to a poor fit. The mixtures of chemical species which 
we have found lead to the best overall fit to the PAO's 
spectrum, (Xmax), and RMS(Xmax) measurements (80% 
nitrogen and 20% iron, 80% silicon and 20% nitrogen, 
and 90% silicon and 10% protons) accommodate the data 
only about as well as in the case of a pure silicon injection 
spectrum. Thus, the consideration of a mixed, two com- 
ponent, chemical composition from cosmic ray sources 
does not significantly improve the agreement with the 
PAO data. 



IV. THE IMPACT OF INTERGALACTIC 
MAGNETIC FIELDS 

A charged cosmic ray moving through a uniform mag- 
netic field will be deflected by an angle 9o = LcoIi/Rl, 
where L^oh is the size of the region being traveled through 
(the coherence length) and Rl is the Larmor radius of 
the particle. Over the course of its propagation, a typi- 
cal UHECR will travel through a large number of regions 
with differing magnetic orientations, which over a dis- 
tance D leads to an overall deflection of « y^D/Lcoh- 
For nano-Gauss scale flelds with coherence lengths of 
1 Mpc, this yields 




where E is the energy of the cosmic ray, Z is its charge, 
and B is the strength of the intergalactic magnetic field. 
Such deflections cause cosmic rays to propagate over 
longer effective distances than the straight line distance, 
with the path extension, in the small angle limit, being 
approximately given by: 
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(5) 



Magnetic deflections impact particles with the greater and composition that are predicted to reach Earth [3]. 
electric charge and/or lower energy more than other cos- In Figs. [5] [6j and [7j we show our results for all-iron, 
mic rays, and thus effect the resulting spectral shape all-silicon, and all-nitrogen accelerating sources, respec- 
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FIG. 4: The same as in previous figures, but for selected cases in wliicli a mixed cliemical composition is injected from cosmic 
ray accelerators. From top-to-bottom, the chemical mixture at injection is given by: 80% protons and 20% iron, 80% nitrogen 
and 20% iron, 30% nitrogen and 70% silicon, 20% protons and 80% nitrogen, and 50% protons and 50% silicon, respectively. 
The injected spectral index, a, was selected in each case to best match the spectral shape reported by the PAO. The ratios of 
different elements were chosen to provide the best fit to the combination of (Xmax) and RMS(Xmax). In each case, we have 
used -Emax of either 10^°'^ eV or lO'^^ eV, depending on which yields the best fit, where -Emax refers to the cutoff energy for iron 
nuclei, which is related to the cutoff for a species with electric charge Z by Em!i^,z = {Z/2&) x -Emax. See text for more details. 
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tively, and for intergalactic magnetic fields of 0.0, 0.1, and 
0.3 nG, with coherent scales of 1 Mpc for the non-zero 
B-field cases. In the case of all-iron injection, the tension 
between the PAO's (Xmax) and RMS(Xmax) measure- 
ments is reduced to some degree by the effects of the mag- 
netic fields, but the spectral fit becomes worse as the field 
strength is increased. In the all-silicon and all-nitrogen 
cases, quite good agreement can be found with the PAO 
data. In particular, for --(B/0.3 nG)x(Lcoh/l Mpc)^/^, 
we find an overall per degree-of-freedom of 1.45 and 
1.14 for the all-silicon and all-nitrogen cases, respectively 
(using EPOS, which yields a better fit than the other 
hadronic simulations). These represent far better fits 
than we have found without magnetic fields. 

There are also mixtures of injected nuclei that can pro- 
vide excellent fits to the PAO's data if the effects of inter- 
galatic magnetic fields are included. In Fig. [8j we show 
results for the two mixtures that we found to provide 
the best overall fits to the PAO's data. As expected, a 
mixture of injected nitrogen and silicon provides a very 
good fit (x^ per degree-of-freedom less than 1.0 for mix- 
tures with between 10% and 50% sihcon). We also find 
that injected nitrogen with a small 10%) admixture 
of iron yields a very good fit with 0.3 nano-Gauss mag- 
netic fields with coherence lengths of 1 Mpc. In each 
case considered here, including an admixture of protons 
in the injected flux from UHECR sources only worsens 
the overall quality of the fit. 

To some extent, the key aspects of these results follow 
from the consideration of the unfolded arriving UHECR 
composition at Earth from the PAO data. As an exam- 
ple case to demonstrate this, we consider the ^lO^^ '^ eV 
(ATmax) and RMS(A,nax) PAO data point. For the EPOS 
model, this point is found to sit on the line expected for 
nitrogen nuclei for both (Xmax) and RMS(Ainax)- The 
possibility that a mixed composition may be also con- 
sistent with the data is then tested through the consid- 
eration of fractional admixtures, in multiples of 1/5, of 
protons, helium, nitrogen, and iron. Only the case of a 
small admixture of 20% helium is found to still be able to 
provide agreement for both the (Amax) and RMS(Xmax) 
data (staying within the error bars of their respective 
data point), no other possible admixtures are found to 
provide such agreement with this data point in both data 
sets. Thus, in order to ensure such a composition arrives, 
the injected composition must be both relatively narrow 
(in charge distribution) and heavier than nitrogen. 



V. DISCUSSION AND CONCLUSIONS 

In this article, we have studied the propagation of 
ultra-high energy cosmic ray nuclei in an attempt to 
find scenarios that provide a good fit to the observa- 
tions of the Pierre Auger Observatory (PAO), including 
the measured spectrum, average depth of shower maxi- 
mum ((Amax)); and the RMS fluctuations of the depth 
of shower maximum (RMS(A',„ax))- 



We have found that for all values of injected maximum 
energy and spectral index considered for the UHECR 
sources, the PAO's data, for energies above 10^^'^ eV, 
is best accommodated in cases in which a homogeneous 
distribution of sources inject largely intermediate mass 
nuclei, ranging from roughly nitrogen to silicon. This 
result is found for all hadronic models considered, with 
the best fit results being found using the EPOS model. 
Furthermore, within this same range of source parame- 
ters, the presence of ~(B/0.3 nG)x(Lcoh/l Mpc)^''^ in- 
tergalactic magnetic fields lead to far better fits than are 
found in the absence of cosmic ray deflection. No good 
fits to the PAO's measurements were found without sig- 
niflcant intergalactic magnetic fields.^ Along with these 
results we note that for the EPOS model, whose results 
are responsible for our best fits, no good fits were found 
with the in any scenario considered in which the sources 
of the UHECRs inject a significant (> 10%) fraction of 
protons. Should the PAO data points all be shifted up 
in energy by approximately 20%, as has been suggested 
[30] . the best fit results, for all cases quoted here, are 
found to shift towards a heavier nuclei composition. 

These results have a number of interesting and impor- 
tant consequences. Firstly, if the cosmic ray spectrum 
consists largely of intermediate mass nuclei at the high- 
est energies, as our results suggest is the case, one would 
not expect any significant correlation between the arrival 
directions and sources of such particles. For a /^G-scale 
Galactic Magnetic Field with a correlation lenth of ^1 
kpc, a 6 X 10^^ eV nucleus with Z « 10 would be deflected 
by on the order of 10°, in addition to any deflection re- 
sulting from intergalactic flelds. This conclusion, per- 
haps, suggesting difficulties for the prospects of UHECR 
astronomy, for all but the most bright local sources. 

Under the assumption that the UHECRs consist 
uniquely of protons, a "guaranteed" flux of UHE neu- 
trinos is predicted to be produced as a result of their 
propagation. This signal, known as the cosmogenic neu- 
trino flux [31j , is potentially within the reach of experi- 
ments such as IceCube ANITA 33J, and the PAO 
itself |33| . If the the UHECR spectrum consists of heavy 
nuclei, however, the cosmogenic flux can be considerably 
suppressed [35]. In the scenarios that we have found to 
fit the PAO's data, the cosmogenic neutrino spectrum is 
predicted to be suppressed by an order of magnitude or 
more relative to the all-proton case [36' . The neutrino 
flux from the sources of the UHECRs themselves is also 
expected to be suppressed in these scenarios [37]. 

The authors note that since the submission of this 
work, new results on the UHECR composition from 
(Amax) measurements have been published from the 
HiRcs experiment [5S] . These results appear at first sight 
to sit in contradiction to those from the PAO. Along with 



^ Although magnetic fields in the vicinity of the sources could also 
plausibly account for this conclusion. 
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FIG. 5: Results for an all-iron injection spectrum with negligible magnetic fields (top), 0.1 nano-Gauss (middle), and 0.3 
nano-Gauss (bottom) intergalactic magnetic fields. For the non-zero B-field cases, a coherence length of the magnetic field of 
1 Mpc has been used. 



the apparent contradiction between these data sets, other 
recent results from the Yakutsk experiment fail to re- 
solve this confusion, with different analysis of the data 
arriving at different results for the UHECR composition 
[321 130], an issue whose origin may lie in the different 
methods used to asign energy to the events gT]. Due to 
the strength of the conclusions that have been possible 
from PAO results, the conflict between these sets of re- 
sults are considered a pressing issue which we hope the 
experimental community will thoroughly investigate. 
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